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ABSTRACT 

o " 

' Context. All Globular Clusters (GCs) studied in detail so far host two or more populations of stars (the multiple population phe- 

nomenon). Theoretical models suggest that the second population is formed from gas polluted by processed material produced by 
massive stars of the first generation. However the nature of the polluter is a matter of strong debate. Several candidates have been 
proposed: massive main-sequence stars (fast rotating or binaries), intermediate-mass AGB stars, or SNell. 

Aims. We studied red giant branch (RGB) stars in the GC M4 (NGC 6121) to measure their chemical signature. Our goal is to mea- 
l/^ ' sure abundances for many key elements (from Li to Eu) in order to give constraints about the polluters responsible for the multiple 

populations. 

Methods. We observed 23 RGB stars below the RGB-bump using the GIRAFFE® VLT2 spectroscopic facility. Spectra cover a wide 
range and allowed us to measure Mght (Li,C,'^C/'^C,N,0,Na,Mg,Al), a (Si,Ca,Ti,) , iron-peak (Cr,Fe,Ni), hght-s (Y), heavy-s (Ba), 
and r (Eu) elements. We completed this study by analyzing a subsample of the UVES spectra presented in Marino et al. (2008) in 
order to have further clues about light s-elements of different atomic number (Y and Zr). 
(~| , Results. We confirm the presence of a bimodal population, first discovered by Marino et al. (2008). Stars can be eas- 

^ |i ily separated according to their N content. The two groups have different C,'-C/''C,N,0,Na content, but share the same 

I ' Li,C-l-N+0,Mg,Al,Si,Ca,Ti,Cr,Fe,Ni,Zr,Ba and Eu abundance. Quite surprisingly the two groups differ also in their Y abundance. 

O ' This result is strongly supported also by the analysis of the UVES spectra. 

Conclusions. The absence of a spread in a-elements, Eu and Ba makes SNell and AGB stars unlikely as polluters. On the other hand, 
^2 ' massive main-sequence stars can explain the bimodality of Y through the weak s-process. This stement is confirmed independently 

^ , ' also by literature data on Rb and Pb. The lack of a Mg/Al spread and the extension of the [0/Na] distribution suggest that the mass of 

the polluters is between 20 and 30 Mg. This implies a formation time scale for the cluster of 10-r30 Myrs. This result is valid for M4. 
Other clusters like NGC 1851, M22, or co Cen have different chemical signatures and may require other kinds of polluter. 

^ , Key words. Galaxy: Globular Cluster:individual: M4 - stars: abundances, light and s-element content 

a^ 

• 1 . Introduction the older generation which have experienced hot H-burning via 
■ p-capture. In this picture the older generation is the most He-N- 

CTN In the last few years, following the discovery of multiple popula- ^a-Al poor and C-O-Mg rich, while the younger generation is 

^ tions m the color-magnitude diagrams (CMD) of some globular ^^^^^^^ an enhancement of its He content, together with N, 

^ clusters (GCs) and m spectroscopic samples of many of them, ^a, and Al, while C, O, and Mg turn out to be depleted. This hy- 

. . the debate on their formation has been renewed. In this respect, p^thesis can also explain correlations or anti-correlations of light 

> the most interesting and pecuhar clusters are w Centauri and elements present at the level of unevolved stai's (.Gratton et all 

• NGC 2808, w here at least 3 main sequences (MS) are present ifooH) 
rS ' (iBedin et al. f2004: Villanova et al. 2007; Piotto et al. 2007). 

h If those two objects represent the most extreme cases. Pollution must come from more massive stars. The main 

- - ' it is now recognize d that all GCs studied in detail so far classes of candidate polluters are: intermediate-mass AGB stars 

(ICarretta etaI.ll25n^ show at least some kind of spread in their ([Ventura et al. || 2004 4<M<7 M^), fast-rotating massive main- 



CD 



light element content at the level of the RGB, the most ev- sequence (MS) stars (IPecressin et al. 1 12007. M>15 Mo), and 

ident bein g the spread in Na and O, elements that are anti- ^^^o massive MS binary stars (de Min k et al. || 2009t M~20 M^). 

correlated (ICarretta et al. lllml . Na and O abundances are also ^11 these channels can potentially pollute the existing interstellar 

(an ti)correlated w ith other light elements, such as C,N,Mg, and material with products of complete CNO cycle where N is pro- 

AI jGratton et al 1120041) duced at the expense of C and O, the NeNa cycle, where Na is 

The most natural explanation for this phenomenon is the Pi^duced at the expense of Ne, and also the MgAl cycle where 

self-pollution scenai-io, where a cluster experiences an extended produced at the expense of Mg (see Renzini 2008 for an 

cx.tcrisivc review I 

period of star formation, with the younger population born from ^' 

an interstellar medium polluted by ejecta coming from stars of AGBs eject part of their outer envelope during the thermal 

pulses after they undergo hot-bottom burning, while massive 

* Based on observations made with ESO telescopes at La Silla stars eject material through stellar winds contaminated by pro- 

Paranal Observatory under program ID 083.B-0083 cessed material brought to the surface because of the fast rota- 
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Fig. 1. The CMD of M4 with the observed RGB stars indicated 
as filled circles. 



produce both light ( i.e. Rb, Sr, Y, Zr) a nd heavy (i.e. Ba, La , 
Ce, Nd) s-elements (iBusso et al. 1 120011: iTravadio et al. [[2004 
through the main-s process, while massive MS stars (M>15M0) 
produce only light s-elemen t (up to A~90, i.e. up to Y or Zr) 
through the weak-s process dRaiteri et al. 11199^ . Finally SNell 
produce mainly a-elements (e.g. Si and Ca) as well as r elements 
(e.g. Eu), while SNela produce mainly iron-peak elements (e.g. 
Fe and Ni) beside r elements ( Wallerstein et al. 1997)). We can 
see that the study of a, Fe-peak, Ught and heavy s, and r elements 
is crucial to disentangle the proposed scenarios. 

The aim of this paper is to measure abundances for a large 
sample of elements in RGB stars of M4 (NGC 6121) in order 
to help constrain the n ature of the pollu t ers. This cluster has 
been studied in d etail divans et al. 1 119991: iMarino et aO 120081: 
lYong et al.1 l2008'). It has a bimodal Na-O distribution. For this 
reason and because of its proximity, it is the ideal target for our 
purposes. We would like to verify which one of the proposed 
polluters, if any, can better explain the observed abundance pat- 
terns. We will focus on a variety of s-elements, but we will in- 
clude in our analysis also many other elements (from Li up to 
Eu) in order to chemically characterize the two sub-populations. 

In Section 2 we describe the observations. In Sec. 3 and 4 
we discuss the determination of the abundances and present the 
results. In Sec. 5 we discuss the results, while Sec. 6 gives the 
conclusions. 



tion or binary interaction. In both cases the primordial material 
(and the older generation) has the same composition as Galactic 
Halo field stars (i.e. He-N-Na-Al poor and C-O-Mg rich), while 
the contaminated material (and the younger generation) is He-N- 
Na-Al rich and C-O-Mg poor. The material required to form the 
second gener ation is kept in the clu ster due to the strong gravi- 
tational field dP'Ercole et al. II2008I) . 

Another scenario was proposed by iMarcolini et aL I (|20Q9|). 
According to this paper in a primordial metal-poor medium pre- 
polluted by SNe II explosion, AGB stars start to eject their en- 
velopes and a simultaneous SN la explosion further contami- 
nates (mainly with iron-peak elements) and collects the ejecta in 
a central region. Here a first (older) generation is formed that, 
at odds with the models described before, is He-N-Na-Al rich 
and C-O-Mg poor. After that the most massive stars of this first 
generation evolve and explode as SNell, that pollute the remain- 
ing gas (mainly with a-elements) and mix it with the primor- 
dial medium. From this new material a second (younger, He- 
N-Na-Al poor and C-O-Mg rich) generation is formed. SNell 
produce also some amount of iron-peak elements but, due to this 
mix with the primordial metal-poor medium, stars of the sec- 
ond generation have the same iron con tent as the first genera- 
tion. For our purposes the main point of IMarcolini et al. rd2009l) 
is that SNell are responsible for the chemical inhomogeneities 
observed nowadays. 

Also the abundance of other elements (including s- and r- 
process elements) may differ in stars of the first or second gen- 
eration according to the nature of the polluters. 

While the pollution scenario is widely accepted nowadays, 
a key piece of information is still missing. We need to verify 
which kind(s) of polluter is responsible for the contamination. 
This will help constrain the conditions and timescale of the pol- 
lution process. 

All processes described above produce (anti)correlations in 
light elements (from C up to Al), but they behave differently as 
far as other elements are concerned. AGB stars are known to 



2. Observations and data reduction 

Our dataset consists of high resolution spectra collected in June- 
August 2010. The spectra come from 10x45m exposures, ob- 
tained with the FLAMES -GIRAFFE spectrograph, mounted at 
the VLT telescope. Weather conditions were good with a typical 
seeing of ~1.0 arcsec. We selected 23 isolated stars at V~ 14.5, 
loc ated below the RGB-b ump of the cluster, from B, V photome- 
try dMomanv et alj2003l, see Fig.[T]). All the stars He within +0.3 
mag, so can be assumed to be in the same evolutionary phase. All 
stars were observed with 4 different set-ups, HR04 (4 exposures, 
ranges 4188-4392 A,R=20000), HRll (2 exposures, range= 
5597-5840 A, R=24000), HR13 (2 exposures, range =6 120-6405 
A, R=22000), and HR15N (2 exposures, range= 6470-6790 A, 
R= 17000). 

Data were reduced using the dedicated pipeline 
BLDRS vO.5.3, written at the Geneva Observatory (see 
http://girbldrs.sourceforge.net). Data reduction includes bias 
subtraction, flat-field correction, wavelength calibration, sky 
subtraction, and spectral rectification. Spectra have a typical 
S/N of -150 at 6300 A. 

Radial velocities were measured by the fxcor package in 
IRAF, using a synthetic spectrum as a template. The mean helio- 
centric value for our t argets is 71.9±0.9 km/s, while the disper- 
sion is 4.2+0.6 km/s. ISommariva et al. I d2009l) gives 70.3+0.2 
km/s as heliocentric radial velocity for M4, a nd the typical dis- 
persion for a cluster of its mass is ~4-5 km/s dPrvor & Mevlan I 
1993), in very good agreement with our results. On the basis of 
this result, we conclude that all our targets are cluster members. 

Table [1] lists the basic parameters of the selected stars: the 
ID, the J2000.0 coordinates (RA & DEC), U,B,V magnitudes 
(Momany, private communication), heliocentric radial velocity 
(RVh), Teff, log(g), micro-turbulence velocity (vt). For determi- 
nation of atmospheric parameters see the next section. 
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Fig. 2. Example of spectral synthesis for Y and Ba lines applied to GIRAFFE data (left panels, star #907), and to UVES data (right 
panel, star #19925). Abundances used in the spectral synthesis are indicated. 



Table 1. Basic parameters of the observed stars. 



ID 


RA(degrees) 


DEC(degrees) 


U(mag) 


B(mag) 


V(mag) 


RVH(km/s) 


Teff(K) 


log(g)(dex) 


v,(km/s) 


28590 


245.81437500 


-26.64647222 


16.13 


15.60 


14.52 


73.8 


4920 


2.65 


1.07 


33584 


245.79137500 


-26.47600000 


16.14 


15.60 


14.50 


77.0 


4940 


2.67 


1.07 


36820 


245.98212500 


-26.65019444 


16.33 


15.81 


14.69 


71.0 


4940 


2.73 


1.12 


37614 


245.93579167 


-26.63430556 


16.13 


15.68 


14.55 


73.9 


4940 


2.55 


1.10 


39100 


245.95291667 


-26.60716667 


16.04 


15.54 


14.40 


70.7 


4890 


2.68 


1.00 


40197 


245.93508333 


-26.59366667 


16.33 


15.78 


14.65 


79.6 


4940 


2.75 


1.16 


41863 


245.86783333 


-26.57441667 


16.40 


15.83 


14.71 


67.4 


4940 


2.65 


1.24 


42561 


245.95445833 


-26.56680556 


16.22 


15.70 


14.59 


66.4 


4950 


2.77 


0.94 


43020 


245.87895833 


-26.56230556 


16.48 


15.99 


14.87 


76.1 


5030 


3.05 


1.08 


43085 


245.86550000 


-26.56161111 


16.33 


15.74 


14.59 


76.4 


5000 


2.95 


1.13 


43494 


245.87850000 


-26.55780556 


16.28 


15.79 


14.66 


65.8 


4970 


2.67 


1.20 


43663 


245.84891667 


-26.55608333 


16.23 


15.68 


14.58 


78.0 


4960 


2.90 


1.00 


45171 


245.86562500 


-26.54222222 


16.26 


15.71 


14.58 


77.6 


4930 


2.90 


0.88 


45200 


245.92258333 


-26.54197222 


16.36 


15.78 


14.65 


67.4 


5010 


2.70 


1.17 


46201 


245.93333333 


-26.53383333 


15.98 


15.47 


14.38 


70.6 


4930 


2.45 


1.08 


47596 


245.86287500 


-26.52325000 


16.07 


15.58 


14.46 


68.4 


4960 


2.70 


1.23 


48499 


245.95400000 


-26.51641667 


15.92 


15.43 


14.31 


76.0 


4960 


2.85 


1.10 


49381 


245.88962500 


-26.50972222 


16.26 


15.78 


14.69 


69.8 


4940 


2.55 


1.22 


50032 


245.93904167 


-26.50455556 


16.13 


15.59 


14.50 


74.0 


5050 


3.03 


0.99 


53602 


245.83920833 


-26.47597222 


16.17 


15.54 


14.42 


68.0 


5000 


3.05 


1.00 


67553 


246.00250000 


-26.47241667 


15.77 


15.26 


14.17 


66.6 


4900 


2.60 


1.02 


8460 


245.94720833 


-26.37500000 


16.14 


15.47 


14.37 


70.4 


4930 


2.50 


1.18 


907 


246.03741667 


-26.37791667 


15.81 


15.31 


14.21 


70.1 


4920 


2.65 


1.13 



3. Abundance analysis 

The chemical abundances for Na, Mg, Al, Si, Ca, Ti, Cr, Fe, 
and Ni were obtained from the equivalent widths (EWs) of the 
spectral lines. See iMarino et al. I (l2008l) for a more detailed ex- 
planation of the method we used to measure the EWs. For the 
other elements (Li, C, N, O, '^C/'^'C, Y, Ba, Eu), whose Hnes ai-e 
affected by blending, we used the spectrum-synthesis method. 
For this purpose we calculated 5 synthetic spectra having dif- 
ferent abundances for the element, and estimated the best-fitting 
value as the one that minimize the rm.s. Na and Al present few 



features in the spectrum, so in this case abundances derived from 
the EWs were cross-checked with the spectral synthesis method 
in order to obtain more accurate measurements. Only lines not 
contaminated by telluric lines were used. 



Atmospheric parameters were obtained in the following way. 
Firs t of all Teg was deriv ed from the B-V color using the relation 
bv Alonso etal.l (Il999l) and the reddening (E(B-V)=0.36) from 
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Fig. 3. Correl ation between light-el ements (C,N,0,Na) as obtained from our GIRAFFE stars. Green points in the upper left panel are 
the results bv lMarino et al. I (12008 *). N-poor stars are indicated as open circles, while N-rich as filled circles. Red crosses represent 
the mean value for each sub-population. 



[Harris I (11996*). Surface gravities (log(^)) were obtained from the 
canonical equation: 

log(^) = log(^) + 4 • log(^) - log(-^) 

where the mass M/Mq was assumed to be 0.8 M©, and 
the luminosity L/Lq was obtained from the absolute magnitude 
My assu ming an apparent distance modulus of (m - M)v = 12.82 
(lHarrisiri996 ). The bolometric co rrection (BC) was de rived by 
adopting the relation BC-Teff from lAlonso et al. 1 (ll999h .Finallv. 
micro-turbulence velocity (vt) was obtained from the relation of 
[Marino etal. (2008). 

These atmospheric parameters were considered as initial esti- 
mates and were refined during the abundance analysis. As a 
first step atmo spheric models were calculated using ATLAS9 
(lKuruczlll970h and assuming the initial estimate of TgfF, \ogig), 
and Vt, and the [Fe/H] value from Harris ( 1996). 
Then Teff, Vt, and \og{g) were adjusted and new atmospheric 
models calculated in an interactive way in order to remove trends 



in Excitation Potential (E.R) and equivalent widths vs. abun- 
dance for Tetf and Vt respectively, and to satisfy the ionization 
equilibrium for log(g). Fel and Fell were used for this purpose. 
The [Fe/H] value of the model was changed at each iteration 
according to the output of the abundance analysis. T he Local 
Thermodynamic Equilibrium (LTE) program MOOG dSneden I 
1973i) was used for the abundance analysis. 



We checked the reliability of our atmospheric parameters by 
comparing photometric and spectroscopic Teff, finding a mean 
difference in temperature lower than 50 K. 
A further check was performed on our log(g) scale. We inverted 
the previous equation in order to obtain the mass and calculated 
the mean mass of our targets. We got M- 0.83 + 0.06 Mq, in 
good agreement with the value obtained from isochrone fitting 
(~0.8 M0). Our conclusion is that our Teff and log(g) values can 
be safely used to obtain abundances. 
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( Gratton et al. ' (2003), 
[McWilHam & Rich ( 1994), McWilliam (1998), SPECTRU 
and SCAIsEI), and calibrated using th e Solar-inverse techni que 
by the spectral synthesis method (see IVillanova et al. Il2009l for 
more details). For this purpo se we used the high resolution, high 
S/N NOAO Solar spectrum (Kurucze t 317111984 '). Adopted solar 
abundances we obtained with our li nelist are reported in Tab. [3] 
and |4] together with those ones by iGrevesse & SauvaTI (Il998l) 
for comparison. We emphasize the fact that all the linelists were 
calibrated on the Sun, including those used for the spectral syn- 
thesis. 

Lines treated with the EQW method are reported in Tab. |6] 
(electronic edition) together with the adopted parameters and 
equivalent widths star-by-star. Parameters for lines treated with 
the spectral synthesis method are not reported because the line 
list would be too long (thousands of lines in some cases). For 
these lines references are given above. 

Li was measured from the line at 6707 A, while our determi- 
nations of C,N,0 abundances are based on the G-band at 4310 
A, the CN band at 4215 A, and the forbidden O line at 6300 A 
respectively. CN lines at 4230 A were used also to estimate the 
i^C/i^C ratio. 

These features were also check ed on the high resolut ion, 
high S/N spectrum of Arcturus. See IVillanova et alTI (1201 Ol) for 
more details. Abundances for C, N, and O were determined all 
together in an interactive way in order to take into account any 
possible molecular coupling of these three elements. 

Our targets are objects evolved off the main sequence, so 
some evolutionary mixing is expected. This can affect the pri- 
mordial C,N,0 abundances separately, but not the total C-hN+O 
content because these elements are transformed one into the 
other during the CNO cycle. Is does not affect the relative C,N,0 
of our stars either, since the stars are all in the same evolutionary 
phase. 

Na was obtained from lines at 6154 and 6160 A (EQW) 
and 5682 and 5688 A, and corrected for NL TE effects follow- 
ing the prescription bv iGratton et al. I (Il999l) . Mg was obtained 
from the line at 5711 A while Al from the lines at 6696 and 
6699 A. Y abundance was measured using the line at 4375 A 
while for Ba we used the line at 6494 A(see Fig. |2] left pan- 
els). For Ba we took the h yperfine splitting into account using 
iMcWiUiam & Rich I ( 11994 andlMcWilliam1( ll998h data. Finally 
Eu was obtained from the line at 6645 A. 

In order to extend and confirm our results, w e analyzed also 
a subsample of 24 stars of iMarino et al. 1 (I2008h observed with 
UVES, the high-resolution spectrograph mounted at VLT tele- 
scope. As parameters we used those published there, but we ex- 
tended the chemical analysis to Y and Zr (see Fig.|2] right pan- 
els), two elements not considered in that paper. Y was obtained 
from the line at 4900 A, while Zr from the line at 6127 A. Also 
in this case abundances were obtained by spectrum-synthesis. 
In UVES spectra two other Y lines were available, at 4883 and 
5087 A, but the one at 4900 A turned out to be stronger and less 
affected by noise. It is partially blended with another line as is 
visible in Fig.|2l but this can be easily managed by the spectrum- 
synthesis method we applied. 



and 



' See 'http://vald.astro.univie.ac.at/~vaId/php/vald.phpl 
http://physics.nist.gov/PhysRefData/ASD/lines_form.html 

See http://w ww.phys. appstate.edu/spectrum/spec trum.html, and ref- 
erences therein 

' See http :// www. astro.ku. dk/~uffea7| 



An internal eiTor analysis was performed by varing Tefj, 
log(g), [Fe/H], and Vt and redetermining abundances of star 
#33584, assumed to represent the entire sample. Parameters 
were varied by ATeff=+50 K, Alog(g)=-i-0.10, A[Fe/H]=+0.05 
dex, and Avt=+0.1 km/s. This estimation of the internal errors 
for a tmospheric parameters was performed as in iMarino et al. I 
('2008). Results are shown in Tab. |5] including the error due to 
the noise of the spectra. This error was obtained, for elements 
whose abundance was obtained by EQWs, as the average value 
of the errors on the mean as given by MOOG, and for elements 
whose abundance was obtained by spectrum-synthesis, as the er- 
ror given by the fitting procedure. cr,o, is the squared sum of the 
single eiTors, while o-Qbs is the mean observed dispersion of the 
two sub-populations of the clusters (as identified by their N con- 
tent, see next section). The agreement between the two values is 
reasonable for all elements. 



4. Results 

First of all, in Fig. [3] we plot the abundance of C,N,0,Na. In all 
the panels we clearly see that the distribution is bimodal for all 
the f our elements consider ed. Na-O anti-correlation is compared 
with IMarino et al.) (l200 8l) (green points). We confirm that stars 
in the cluster are divided in two well separated groups having 
different light-element content. In the following analysis we take 
as reference the lower left panel of Fig.[3]because N appears to 
be the best element to separate the two groups. From now on 
in all the figures and in the discussion of the results, we divide 
our stars into N-poor (open circles in all the figures), and N-rich 
(filled circles in all the figures). 

We note that in our sample 9 stars belong to the N-poor group 
(the so called first generation), while 14 to the N-rich group (the 
so called second generation). This means that about 40% of our 
stars belong to t he first generation. T his a gree within the error s 
(a- ~ +6%) with lCarretta et al. I(l2009h and lMarino et al. I(l2008h . 
where the authors find that the cluster is composed of ~30% and 
~50%, respectively, first generation stars. 

Errors on the single measurements (the black errorbars in 
the figures) for a given element are the (Xobs of Tab. |5] The red 
crosses in the figures represent the mean abundance and the eiTor 
of the mean for each group. We see that N-poor stars are also C- 
rich, O-rich, and Na-poor, while N-rich stars are also C-poor, 
O-poor, and Na-rich, in accord with the theoretical expectations 
outlined previously. 

Mean abundances we obtained for the two groups and for the 
cluster are summarized in Tab.|2] while abundances for each star 
are summarized in Tab.[3l and|4l 

With respect to the mean abundances of the two groups, for 
each element in the 4''' colum of Tab. |2] we report the abundance 
difference (i.e. the significance) in units of (Tei that is defined as: 



el,N- 



-poor el,N- 



rich 



where o-gi^N^poor and (Tei^N-rich are the errors on the mean 
abundance of the two groups as given by the 2"'' and y' columns 
of Tab. |2] This value tell us if this difference is significant with 
a value of cr^/ >3 implying strong significance. The second part 
of Tab. |4] reports the results ob tained from the analysis of the 
UVES spectra of IMarino et al. I ((2008). Na abundances are ob- 
tained from that paper For these stars we do not have the N 
content. However it is clear from Fig. |3] that N-poor and N- 
rich stars can be easily identified also by their Na content. So 
UVES stars were classified according to their [Na/Fe] value. All 
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Table 2. Mean abundances of the two groups of stars (2"'' and 3'''' column). The 4''' column is the the significance of the abun- 
dance difference of the two groups, in units if cTp/. 5''' column gives the mean abundance of the cluster as the average of the two 
group s. 6''\ l"\ and 8'^' col umns are the abundances by Marino et al. (2008, Ma08), Iva ns et al. (1999. Iv99), and Yong et al] 
(l2008l) 4 fYong et alH (l2008bl) =Yo08. respectively. The results of the last two papers are reported together because they are comple- 
mentary. 
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Fig. 4. [N/Fe]->2Qi3(- anti-correlation for our GIRAFFE stars. 
Red crosses represent the mean value for each sub-population. 



stars with [Na/Fe] <0. 23 were considered N-poor, all stars with 
[Na/Fe]>0.23 were considered N-rich. 

Comparing values in Tab.|2]we see immediately that the two 
sub-populations have the same content (difference of 1.8cr in 
the worst case) of a (Mg,Si,Ca,Ti) and iron-peak (Cr,Fe,Ni) el- 



ements. Also the Li and Al abundances and the total C+N+0 
content are the same within the errors. 

On the other hand the light elements C,N,0,Na are different 
between the two groups, with a significance of more than 4 cr. 
Also the carbon isotopic ratio is different. This is shown in Fig.|4] 
were we plot [N/Fe] vs. '^C/'^'C. In this figure an anti-correlation 
appears, and the difference in '^C/'^C between the two groups is 
more than 3 cr, although there is a substantial overlap of the two 
distributions due in part to the measurement error. This result is 
not unexpected because the N-rich population is supposed to be 
born from material more chemically evolved with respect to the 
N-poor one. So it should have a lower carbon isotopic ratio, as 
we find. 

Fig. |5] (upper panels) displays the abundances of the s- 
elements Y and Ba for our GIRAFFE data. As in the case of the a 
and iron-peak elements, the two sub-populations have the same 
Ba content. However the mean Y abundance is significantly dif- 
ferent, at a level of more than 3 cr (see Tab. |2]). We fitted to the 
points a straight line and calculated the slope and its error. Values 
are reported inside each panel. The Ba slope is compatible with 
(i.e. the same mean [Ba/Fe] value for the two groups), while the 
Y slope is incompatible with with a confidence of more than 
3 cr (i.e. two different mean [Y/Fe] values for the two groups, 
although there is substantial overlap between the two distribu- 
tions due in part to the measurement error). We cannot rule out 
a trend, but we favor a bimodality. 

Lower panels display the abundances of the s-elements Y 
and Zr for the UVES data vs. [Na/Fe]. Again the mean Y con- 
tent of Na-poor (N-poor) and Na-rich (N-rich) stars is very dif- 
ferent, at a level of more than 3 cr, while they share the same Zr 
abundance. 

The Y abundance from GIRAFFE and UVES obser- 
vations deserves a further comment. For N-rich stars the 
two databases give very good agreement (+0.31 vs +0.33 
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Fig. 5. Upper panels: [Y/Fe] and [Ba/Fe] as a function of [N/Fe] for our GIRAFFE stars. Lower panels: [Y/Fe] and Zr/Fe as a 
function of [Na/Fe] for out UVES stars. A straight line was fit to the data of each panel and the slope with its error is reported. Red 
crosses represent the mean value for each sub-population. 



dex), well within Icr. Values for N-poor stars instead ap- 
pear different ([Y/Fe]=-H0.10+0.06 dex for GIRAFFE and 
[Y/Fe] =-hO. 19+0.03 dex for UVES). However if we consider the 
errors, the diff'erence of 0.09 dex is significant at the level of 1.3 
cr, too low to imply a real diff'erence, so we can safely attribute 
it to measurement errors. 

The total error on Y (0. 12 dex, see Tab.|5]l due to atmospheric 
parameters and S/N is high. It is dominated by errors in gravity 
and microturbulence, but also S/N gives a not-negligible contri- 
bution. However this is a random internal error, and for this rea- 
son it is fully included in the error of the mean Y abundance of 
the two groups of stars that we use to calculate the significance 
in Tab. HI 

Finally in Fig.|6]we report [Eu/Fe] as a function of [N/Fe]. At 
odds with the Y-abundance discussed before, Eu does not show 
any trend, confirming the result reported in Tab. |2] 

Thus, we find no significant difference between the 
two groups defined by their N (or Na) abundance, in 
Li,C+N+0,Mg,Al,Si,Ca,Ti,Cr,Fe,Ni,Zr,Ba or Eu, but find very 
significant differences in C,N,0,'-C/'^C,Na, and Y. 



4. 1 . Comparison with iiterature 

Firstly we compare our result with iMarino et al. ] (l2008h (see 
Fig. O. Also in that case the authors find a bimodal Na- 
O anti-correlation. They define a Na-poor (green points with 
[Na/Fe] <0.25) and a Na-rich (green points with [Na/Fe] >0.25) 
population, which correspond to our N-poor and N-rich respec- 
tively. Their Na-rich population have the same mean Na con- 
tent as our N-rich stars, while their Na-poor stars have a mean 
Na content that is slightly higher than that of our N-poor pop- 
ulation. This could be a residual of the NLTE correction. Here 
we show that the bimodality is extended also to C, N,'^C/'^C, 
and Y. As lMarino et al. I (l2008l) . we do not find bimodality in a 
(Mg,Si ,Ca,Ti) or iron -peak elements(Cr,Fe,Ni), nor Ba. In par- 
ticular iMarinoetal.l (|2008) Na -poor and Na-rich populations 
have the same [Ba/Fe] within 1 cr, confirming our finding. Thus 
the two populations have the same abundance as far as these el- 
ements are concerned. 

Another important paper is llvans et aTl ( Il999l) . These au- 
thors include also N and indirectly C in their results. They find 
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the same Na-O anticorrelation as we do (see their Fig. 14, upper 
panel). They find also a relatively well-defined C-O correlation 
and N-O anticorrelation as is implied also by our Fig. [3] Their 
C+N+O content is constant within the errors for all the stars, 
but a bit higher that our result (8.24 vs. 8.16 respectively) in ab- 
solute value. As in our case a bimodality is suggested by their 
measurement of the strength of the CN band at 7874 A (see their 
Fig. 1 1), and partially visible in the N-0 anticorrelation where a 
discontinuity at loge(N)~7.70-7.80is visible. 

As suggested by the referee, we comment in more detail on 
the Al content. Marino et al. (2008) found a possible spread in 
Al, confirmed by I vans et al. I (1199 9). where Al is correlated with 
Na with a significance of more than 3cr. We instead find that the 
Al content is the same for the two groups of stars, which implies 
no Al-Na correlation. A possible explanation could be the pres- 
ence of some unrecognized molecular line (CN?) that is blended 
with Al lines . This Une could pre fere ntially affbct colder stars, 
and both the Marino et aTl (l2008l) and llvans etal.l (11999) stars 
are colder than ours. Because C and N have different abundances 
for the two groups, we expect that this hypothetical line has a 
different strength (assuming the same Te/f) if a star belongs to 
one group or the other. This would cause a spurious correlation 
between Al and N or Na also if the Al content is the same for 
the two groups. This hypothesis is not totally unreasonable be- 
cause the spectral ATLAS of ArcturusQ shows that the Al lines 
are possibly blended with some CN line (i.e. the CN line at A- 
6698.746 A), but the apparent lack of an increasing trend in Al 
abundance with evolutionary phase in previous works could ar- 
gue against it. 

Another explana tion could be the e volutionary state of 
iMarino et aD (12008^ and 'Ivans et all (11999 ) targets. Their stars 
are all above the RGB-bump and so aff'ected by more evolu- 
tionary mixi ng. This is also proved by the '^C/'^'C value ob- 
tained by Iv ans et al. I (|l999) that is of the order of 4-5. Our stars 
instead have ^^C/^^C~20, implying a less dramatic mixing. A 
deeper mixing after the RGB-bump could alter significatively 
(and differentially, with the N-rich group being the most af- 
fected) the primordial Al content through products of the MgAl- 
cycle. Howver if nucleosynthesis products of the MgAl-cycle 
are indeed mixed up at this evolutionary phase, this would have 
consequences also for the Na abundances, as the NeNa-cycle op- 
erates at lower temperatures dCharbo nnel 2005 |), and there are 
studies that argue against such changes (Gratton et al. 2000). A 
further discussion of this point is beyond the scope of this paper. 
In any case, in the present paper we find that N-poor and N-rich 
stars have the same Al content. 

Our main results are independent of the absolute abun- 
dances, but we want to add further comments on this 
point. Comp a rison with absolute abun d ances published b y 
IMarino et all (l2()08h. Ilyans et al. I (Il999h . lYong et al. I (l2008h . 
and lYong et al. J (|2008b')~ (the four most recent papers) are re- 
ported in Tab.|2] The last two papers are reported as one (Yo08) 
in the table and in the following discussion because they are 
complementary. For some elements the mean abundance of the 
cluster agrees well (difference of 0. 1 dex or less) with these four 
papers. This is true for C+N+O, Mg, Ti, Cr, and Fe. ^^C/"C 
is higher in our case, as expected by the evolutionary stage 
of our stars. Our Al , Si, and Ba content agrees well also with 
IMarino et al. ] (l2008h . but in this case the scatt er is larger (be- 
tween 0.1 to 0.3 dex) with respect to llvans et al. (19 99h and 
Yo08 Our rO/Fe1, [Na / Fel, a nd [Ni/Fe] match IMarinoet al. I 
(120081) and llvans etal.l (Il999l) within 0.1 dex, while the dis- 



agreement is worse (0.22, ~0.2, and 0.13 dex respectivelly) with 
respect to Yo08. On the other hand our [Ca/Fe] matches Yo08, 
but it is 0.13 and .15 dex higher then Marino et al. (2008) and 
llvans etal.l(ll999l) . The difference in C and Eu is a b it high (0.2 2 
and 0.15 dex respectively) with respect to Ivans et al. I d 19991) . 
but it is large (0.37 dex) only in the case of N. Finally Yo08 ob- 
tained a Zr value 0.23 dex lower than our, a Eu value 0.20 dex 
higher, while the disagreement is large as far as [Y/Fe] is con- 
cerned (-0.4+0.5 dex). 

We further comment about Ba and Eu. W e find 
[Ba/Fe]~+0.3, in agreem ent with IMarino et aTl (|2008|) 
who give [Ba/Fe]~+0.4. Ivans et al. (1999^ gives a much 
higher value ([Ba/Fe]~+0.6) while Gratton et all (Il986l) gives 
[Ba/Fe] — hO.O. Our value is in the middl e of the literature 
values. We find [Ba/Eu]=+0.10, lower than llvans et"an (Il999l) 
but still higher than the solar system value and much higher 
then field halo and globular cluster giants, where [Ba/Eul is 
typica lly negative with a range from -0.2 to -0.6 (llvans et al. I 
Il999h . This result confirms that M4 has a larger s- to r-process 
contribution than in the Sun, and also supports the llvans et aTl 
(1999) suggestion that the period of star formation and mass 
loss that preceded the formation of the observed stars in M4 
was long enough for AGB stars to contribute their ejecta into 
the primordial ISM of the cluster. 
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Fig. 6. [Eu/Fe] vs. [N/Fe] for our GIRAFFE stars. Red crosses 
represent the mean value for each sub-population. 



5. Discussion 

As discussed in the introduction, all GCs studied in detail to date 
show some kind of spread in their light-element (from C to Al) 
abundance. The amount of the sp read varies a lot from cluster to 
cluster Some GCs, such as M22 (Marino et al. l2009h . show also 
a spread in a and iron-peak elements, but this is an uncommon 
feature. So in this discussion we consider only M4-like objects, 
i.e. those clusters having only a spread in light elements (and 
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possibly in s-elements). However M4 appears to be rather unique 
in that the 'spread' is actually a bimodality. 

The most natural explanation for this phenomenon is the 
self-pollution scenario, where a first generation of stars is formed 
from primordial material. In the most accepted model, this mate- 
rial is O-rich and Na-poor with respect to the second generation 
that will form later Then some class of stars (massive MS stars, 
either fast rotating or binaries, or intermediate-mass AGB stars) 
of this first generation pollute the interstellar material. This ma- 
terial (O-poor and Na-rich) is kept in the cluster due to the strong 
gravitational field of the massive cluster, and it gives rise to a 
new generation of stars. In this model the O-rich/Na-poor stars 
are the oldest. Also the abundance of other elements (including 
He or other light and s-process elements) may differ in stars of 
the first and second generation. 

Another scenario postulates instead that the primordial ma- 
terial the first generation will form from is O-poor and Na-rich 
because of initial pollution by SNela and AGB stars. Then SNell 
of the first generation explode polluting the residual material, 
and a second generation is formed, being O-rich and Na-poor. 
The final result is the same, but in this case the O-poor/Na-rich 
stars are the oldest. 

Each one of these polluters has its own chemical signature 
as discussed in the introduction. SNell ejects material that is 
strongly a-enhanced but only slightly enhanced in iron-peak ele- 
ments. The Marcolini et al. (2009) scenario manages to produce 
a second generation that has the same iron-peak element content 
as the first generation because of the mix with primordial Fe- 
poor material. Nothing is said in that paper about a-elements 
with the exception of O and Mg. However the first generation is 
poor in a-elements (Si,Ca,Ti) because it is formed from material 
polluted by SNela and AGBs (that do not produce these elements 
or only in a negligible amount), while the second generation 
must be rich in a-elements (Si,Ca,Ti) both because of the con- 
tamination by the SNell of the first generation, but also because 
of the mix with primordial material that was pre-contaminated 
by previous SNell. 

The exact estimation of the difference in Si,Ca,Ti content is 
beyon d the pu rpose of this pjiper, but we can give a rough num- 
ber. In iMarcolini et al. I (l2009t Tab. 1) we can see that the pro- 
genitor material of the first generation is enhanced by ~1 dex 
in its Fe content with respect to the primordial one, due to the 
SNIa explosion. So the material the first generation will form 
from is dominated by the products of the SNIa ejecta (plus the 
AGB ejecta, that however do not affect Si,Ca,Ti content). The 
second generation instead is formed from material whose abun- 
dance ratio is dominated by SNell ejecta. So we can estimate 
the difference in Si,Ca,Ti between the two generations compar- 
ing the Si,Ca,Ti content of Galactic stars with [Fe/H]<-1.0 dex, 
whose abundance ratio is dominated by SNell ejecta, and stars 
with [Fe/H]~0.0, whose abundance ratio is dominated by SNela 
ejecta. According to Pompeia et al. ( 2008, Fig. 8), the difference 
of the order of A[Si,Ca,Ti/Fe]~0.3. This is our reference value. 
Because SNell produce also r-element (i.e. Eu), the second gen- 
eration should be also Eu-enhanced. 

AGB stars are the main producers of s-elements, including 
both light-s such as Y, Zr, and heavy-s like Ba, through the main 
component of the s-process. 

Massive MS stars (M>15M0, both as fast rotators or in bi- 
nary systems) produce only light s-elements (e.g. Y) through the 
weak-s process but not heavy s-elements. 

We now compare these predictions with our observational 
results. N-poor and N-rich stars have the same Si,Ca,Ti content 
within a few hundredths of a dex. They also share the same Eu 



abundance. According to these resu lts and comparing theni with 
the Si,Ca,Ti content expected from ' Marcolini et al. I (12009 ^ sce- 
nario, SNell are not viable candidates for the polluters responsi- 
ble. 

The two groups also have the same Ba content, suggesting 
that AGB stars cannot be responsible for the pollution. This 
statement is reinforced if we compare our results with theoretical 
predictions. Considering both GIRAFFE and UVES databases, 
we have a differenc e in [ Y/Fe] for the tw o groups of stars of 0.18 
dex. Accordi ng to Busso et al. I (1200 ll) or the more recent pa- 
per by Karak aset al. I ( 2OI0I) . we would have expected to see an 
equal (if not larger) difference in [Ba/Fe], assuming AGB stars 
as polluter s. And we have the ob servational counterpart of that. 
NGC 1851 (I Villanova et al. Il2010h hosts two distinct populations 
that differ in [Y/Fe] at the level of 0. 1 1 dex. They have a differ- 
ence in [Ba/Fe] that is much larger, 0.41 dex. So for NGC 1851 
we can postulate AGBs as the most probable polluters, while for 
M4 it is very unlikely. 

We are left with massive MS stars. These objects should 
pollute the interstellar material with light s-elements (besides 
C,N,0,Na,Al) produced through the weak s-process. We found 
that the two populations differ in their Y content at a level 
of more than 3 cr, but they have the sa me Zr content. This 
results fits with the theoretical scenario (iRaiterietal. 1 119931: 
'Travagl io et al. Il2004t) that says that the weak s-component is 
responsible for a major contribution to the s-process nuclides up 
to A~90 with a peak for A=80 (^"Kr). The last nucUde of the 
chain is uncertain due to the many theoretical uncertainties, but 
it is not unreasonable (and it is compatible with the models) to 
postulate that, according to our results, the last nuclide produced 
in a significant amount is ^^F, leaving the next nuclide, '°Zr, 
and heavier s-elements, relatively scarce. Unfortunately we can- 
not measure s-elements lighter then Y (i.e. Rb or Sr) in order 
to further investigate our sta tement. Ho wever we can obtain an 
important confirmation from lYong et al. I (12008b). These authors 
estimated Na, Rb, Y, and Pb for a sample of targets, and their re- 
sults are reported in Fig.|7]for those stars with all three elements 
measured. First of all also in this case Na appears to be bimodal, 
with a gap between 0.2<[Na/Fe]<0.4. But, most important, both 
Y and Rb appear to have a positive trend with Na (or N), exactly 
as we have found (but only for Y). This would confirm that light 
s-elements like Y or lighter have a different mean abundance in 
the two populations of M4. The trends, taken separately, have 
a significance of only 1.2-fl.3 cr, as shown in Fig. |7] but if we 
compare the mean Rb content of the two groups (red crosses in 
Fig.|71 lower panel) that are of -hO.37+0.02 and -1-0.42+0.02 re- 
spectively, we have a significance of 1.8 cr. We obtain an even 
stronger confirmation if we consider the two trends together and 
apply the following ab absurdum argument. Let's assume that 
our previous result is wrong and that the two populations have 
the same Rb and Y content. With this hypothesis and using the 
Kolmogorov-Smirnov test, the probability of having a Y vs. Na 
data distribution like that one in Fig.|7]is 17%. For the Rb vs. Na 
data distribution we have a probability of 9%. These two values 
match with the significance at 1.2-H1.3 cr found above. Under 
the previous hypothesis the relations visible in Fig. |7] would be 
due only to measurement errors, so they would be independent 
and we can calculate the probability of having both relations we 
see in Fig.|7]by multiplying the two individual values. The final 
result is only 1.5%. So the probability of the initial hypothesis 
to be true is only 1.5%. This means that lYong et al. I (l2008bl) 
confirms our result of a bimodality in light s-elements with a 
confidence level of 98.5%. Another hint comes from Pb. This is 
a pure s-element, produced only in AGB stars dTravaglio et al. I 



9 



S. Villanova and D. Geisler: Bimodality of light and s-elements in M4 (NGC 6121). 



2004. The mean [P b/Fe] abundances of the two Na groups in 



Yong et al. I (l2008bh are +0.26+0.05 and +0.31+0.02 respectiv 



elly. In this case the significance is 0.8 cr. This further rules out 
AGB stars as candidate polluters, in agreement with our state- 
ment based on Ba. Because of this and because we have double- 
checked with two different spectrographs and two different lines 
the Y-bimodality, we are led to the conclusion that the best can- 
didates for the Y-enrichment for the second generation in M4 is 
the weak s-process, which implies that massive (M>15 Mq) MS 
stars are the best candidates for the self-pollution scenario. 

A further conclusion can be obtained from Mg and Al. These 
elements are processed in the Mg-Al cycle where Al is produced 
at the expense of Mg. The fact that we find the same Mg and Al 
abundance for N-poor and N-rich stars implies that the Mg-Al 
cycle did not activate, at least in the ejected material responsible 
for the pollution. This implies an upper limit for the tempera- 
ture: the material burns at a temperature lower than 50x10^ K 
(iDecressin et al. Il2007) . This implies also that the massive stars 
responsible for the pollution must have been less massive than 60 
Mq. Another hint comes from the difference in the mean [O/Na] 
value of the two groups. In our database it is of the order of 0.7 
dex, while iMarino et al. I (|2008|) give 0.5 dex, so we assume a 
mean value of 0.6 dex. lOecressin et al. ] (I2007D (see their Fig. 
10) give the extension of the [O/Na] distribution as a function of 
the mass of the polluter 0.6 dex implies a mass higher then 20 
Mq but significantly lower than 30 M©. A value between 20 and 
30 Mq seems reasonable. 

Our conclusion is that the best candidate in order to explain 
the abundance spread in M4 are massive MS stars (fast rotators 
or binaries) with masses of the order of 20+30 Mq. Interestingly 
enough, there is another cluster with the same chemic al signature 
as M4 , that is NGC 6397, recently investigated by iLind et al7l 
(1201 1*). In that paper the authors find a bimodality in Na for the 
cluster as in M4, but also a difference in [Y/Fe] between the two 
sub-populations at a level of almost 3cr. As in our case the con- 
tent of the other neutron-captured elements (Zr,Ba,Ce,Nd,Eu), 
all heavier than Y, is the same within ~lcr. They say that the 
[Y/Fe] differecence is 'too small to be convincing' (citation), but 
in the view of our result it could be real. Based on the abun- 
dance pattern, they support the MS massive star scenario too. 
These findings suggest also the hypothesis that a bimodality in 
light-elements (C,N,0,Na) could be a normal feature of the MS 
massive star based self-pollution scenario. 

We underline the fact that our conclusion is valid for M4 
(and possibly NGC 6397). Other clusters show a Mg-Al anti- 
correlation or light an d heavy s-element spread (i.e. NGC 1851, 
IVillanova et al. 1201 Oh . so AGB stars can be involved. 

Finally some clusters like co Centauri or M22 have a spread 
in Ca and Fe, so SNell must also be at work. The picture is 
that each cluster may have its own peculiar process of forma- 
tion within the pollution scenario, and that the polluters can be 
different. 

A peculiarity of M4 is its bimodal distribution, at odds with 
man y other clusters wher e the Na-O anticorrelation is continu- 
ous (ICarretta et al. Il2009l) . This implies that the clusters had a 
first star forming burst followed by a quiescent period, where 
processed material started to flow and mix with the primordial 
gas. During this period no star was formed and the gas had time 
to homogenize, otherwise we would have observed a continu- 
ous anticorrelation as in the other clusters. After that, a second 
burst happened with the formation of the second generation. The 
exact time scale of this process is unknown. However it must 
have been longer than the evolutionary time of very massive stars 
(some Myrs), but shorter than the evolutionary time of AGB stars 



(>40 Mvrs. lVentura & D'Antonal2009h . For the 20+30 Mq star 
polluters postulated above, the evolutionary timescale is of the 
order of -10+30 Myrs. 



0.4 



a =0.11 ±0.09, Signiricance=1.2CT 



o.a 0.4 

[Na/FcJ 



0.6 




0.2 0.4 
[Na/Fe] 

Fig. 7. [Y/F e ] vs. [N a/Fe] and [Rb/Fe] vs. [Na/Fe] for stars from 
lYong et alTI (l2008bh . Linear fits and the slopes with their sig- 
nificance (in units of cr) are reported. For [Rb/Fe] vs. [Na/Fe] 
we plotted the mean abundances (red crosses) of the two groups 
of Na-poor ([Na/Fe]<0.30 dex) and Na-rich ([Na/Fe]>0.30 dex) 
stars. 



6. Conclusions 

In this paper we analyzed a sample of 23 stars belong- 
ing to the RGB of M4 (NGC 6121) and observed with the 
GIRAFFE ©FLAMES spectrograph. Our targets are located be- 
low the RGB bump. We complemented our stu dy by analyz- 
ing a subsample of UVES ©FLAMES spectra of iMarino et al. I 
(l2O08h . We estimated abundances for many key elements: 
Li,C,N,0,i2C/i^C,Na,Mg,Al,Si,Ca,Ti,Cr,Fe,Ni,Zr,Y,Ba,Eu. Our 
aim was to look for some hint in order to help solve the prob- 
lem related to the self-pollution scenario in GCs. According to 
this scenario, a cluster experiences an extended period of star 
formation, where the younger populations were born from an 
interstellar medium polluted by products of the CNO, NeNa, 
and MgAl cycles coming from massive stars of the former gen- 
eration. Postulated polluters include: massive MS stars (M>15 
Mq, fast rotating or binaries), intermediate mass AGB stars 
(4<M<7 Mq), or SNell. Each of these candidates has its own 
unique chemical signature. We confirm the presence of a bi- 
modal population, where the two groups of stars can easily be 
separated by their N content. One group (presumably the old- 
est generation) is N-poor, while the other (the younger genera- 
tion) is N-rich. N-poor and N-rich stars have significantly dif- 
ferent C,N,0,'^C/'''C,Na, and Y content, but share the same 
Li,C+N+0,Mg,Al,Si,Ca,Ti,Cr,Fe,Ni,Zr,Ba, and Eu abundances. 
This rules out SNell because they would produce a second gen- 
eration ff-enhanced with respect to the first one (and Eu en- 
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hanced). AGB stars are also excluded because they should pro- 
duce a Zr and Ba-enhanced (in addition to Y-enhanced) second 
generation through the main s-process. We are left with massive 
MS stars. These objects are able to produce the difference in 
light elements we observe but, most important, they can produce 
also light s-elements like Y through the weak s-process. The Y- 
enhancement of the second generation is the most interesting 
result we found. This scenario is supported also by Yong et aL] 
(|2008b). Based in this paper we found that the two M4 groups 
have different Rb content (that confirms massive MS stars as pol- 
luters) but the same Pb content (that further excludes AGB stars 
as polluters). The lack of Mg/Al enhancement and the extension 
of [O/Na] ratio points toward massive stars with 20<M<30 Mq 
as the most likely polluters. 

Our conclusion is that massive stars in the range 20<M<30 
M0 are responsible for the bimodal population in M4. The time 
scale for the formation of the clusters is ~10h-30 Myrs, with two 
well separated bursts that generated a bimodal population. Other 
cluster have different chemical characteristics, so other kinds of 
polluters and star formations histories are required. 
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Table 3. Individual abundances (light, a, iron-peak elements) of the observed GIRAFFE stars. 



ID 


loge(Li) 


[C/Fe] 


[N/Fe] 


[O/Fe] 


[Na/Fe] 


[Mg/Fe 


] [Al/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Cr/Fe] 


[Fe/H] 


[Ni/Fe] 


'-'c/'Jc 


28590 


1.02 


-0.42 


0.89 


0.26 


0.48 


0.48 


0.66 


0.38 


- 


0.26 


- 


-1, 


.21 


-0.05 


18 


33584 


0.94 


-0.21 


0.35 


0.34 


0.01 


0.54 


0.64 


0.45 


0.41 


0.31 


- 


-1, 


.14 


-0.04 


20 


36820 


0.77 


-0.27 


0.78 


0.34 


0.45 


0.49 


0.64 


0.49 


0.42 


0.32 


0.10 


-1, 


.25 


-0.02 


20 


37614 


1.08 


-0.18 


0.26 


0.43 


-0.08 


- 


0.62 


0.50 


0.36 


0.34 


0.02 


-1, 


.16 


0.04 


22 


39100 


0.97 


-0.19 


0.05 


0.37 


0.03 


0.42 


0.44 


0.47 


0.40 


0.36 


-0.03 


-1, 


.08 


0.00 


25 


40197 


0.97 


-0.27 


0.78 


0.40 


0.35 


0.51 


0.51 


0.51 


0.42 


0.28 


-0.06 


-1, 


.14 


-0.03 


22 


41863 


0.91 


-0.38 


0.87 


0.26 


0.57 


0.51 


0.61 


0.48 


- 


0.23 


0.18 


-1, 


.19 


-0.03 


21 


42561 


0.96 


-0.40 


0.95 


0.28 


0.37 


- 


0.62 


0.39 


0.41 


0.31 


-0.06 


-1, 


.20 


0.04 


13 


43020 


1.05 


-0.19 


0.63 


0.45 


0.23 


0.47 


0.49 


0.51 


0.50 


0.28 


0.01 


-1, 


.12 


-0.04 


22 


43085 


1.01 


-0.32 


0.80 


0.31 


0.40 


0.47 


0.50 


0.42 


0.33 


0.27 


-0.05 


-1, 


.08 


0.02 


19 


43494 


0.93 


-0.18 


0.13 


0.48 


-0.11 


0.49 


0.31 


0.43 


0.44 


0.31 


0.00 


-1, 


.14 


-0.05 


22 


43663 


1.00 


-0.28 


0.69 


0.30 


0.31 


0.36 


0.43 


0.45 


0.33 


0.30 


0.06 


-1, 


.11 


-0.07 


19 


45171 


0.80 


-0.13 


0.06 


0.53 


-0.09 


0.32 


0.44 


0.40 


0.44 


0.38 


-0.09 


-1, 


.11 


0.00 


21 


45200 


0.94 


-0.42 


0.87 


0.14 


0.55 


0.51 


0.57 


0.46 


0.39 


0.36 


0.07 


-1, 


.12 


0.00 


12 


46201 


0.94 


-0.40 


0.81 


0.10 


0.41 


0.54 


0.55 


0.31 


0.44 


0.34 


-0.04 


-1, 


.17 


-0.04 


20 


47596 


1.00 


-0.15 


0.08 


0.48 


-0.03 


0.44 


0.54 


0.35 


0.42 


0.31 


0.04 


-1, 


.16 


-0.02 


22 


48499 


1.17 


-0.15 


0.10 


0.53 


0.03 


0.47 


0.60 


0.54 


0.47 


0.39 


-0.03 


-1, 


.15 


0.04 


25 


49381 


0.95 


-0.27 


0.19 


0.35 


0.02 


0.59 


0.47 


0.38 


0.45 


0.42 


0.13 


-1, 


.19 


0.01 


18 


50032 


1.02 


-0.45 


0.85 


0.11 


0.32 


0.43 


0.50 


0.43 


0.33 


0.33 




-1, 


.04 


0.01 


12 


53602 


1.14 


-0.36 


0.72 


0.23 


0.36 


0.36 


0.39 


0.37 


0.45 


0.33 


-0.11 


-1, 


.03 


-0.01 


18 


67553 


0.90 


-0.30 


0.19 


0.30 


0.14 


0.39 


0.49 


0.35 


0.39 


0.29 


-0.05 


-1, 


.16 


0.04 


20 


8460 


0.99 


-0.36 


0.74 


0.12 


0.40 


0.50 


0.50 


0.37 


0.41 


0.40 




-1, 


.16 


-0.04 


14 


907 


0.81 


-0.48 


0.81 


0.21 


0.42 


0.56 


0.46 


0.37 


0.42 


0.31 


0.03 


-1, 


.20 


-0.03 


14 


Sun:our linelist 






8.49 


7.95 


8.83 


6.32 


7.56 


6.43 


7.61 


6.39 


4.94 


5.63 


7.50 


6.26 




Sun:Grevesse & Sauval (1998) 






8.52 


7.92 


8.83 


6.33 


7.58 


6.47 


7.55 


6.36 


5.02 


5.67 


7.50 


6.25 
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Table 4. Individual a bundances (s and r elements) of the observed GIRAFFE and UVES stars. For UVES also Na abundances from 
iMarino et aTl (l2008h are reported. 



GIRAFFE 



ID 


[Y/Fe] 


[Ba/Fe] 


[Eu/Fe] 


28590 


0.29 


0.33 


0.14 


33584 


- 


0.31 


0.18 


36820 


0.38 


0.31 


0.30 


37614 


0.07 


0.34 


- 


39100 


0.11 


0.25 


0.02 


40197 


0.42 


0.28 


0.10 


41863 


0.23 


0.28 


0.17 


42561 


0.47 


- 


0.17 


43020 


0.28 


0.36 


0.09 


43085 


0.34 


0.37 


0.22 


43494 


-0.09 


- 


0.21 


43663 


0.34 


0.26 


0.26 


45171 


0.37 


0.33 


0.13 


45200 


0.22 




0.13 


46201 


0.07 


0.32 


0.21 


47596 


0.01 


0.28 


0.28 


48499 


0.21 


0.33 


0.30 


49381 


-0.11 


0.26 


0.21 


50032 


0.33 


0.39 


0.19 


53602 


0.42 


0.37 


0.44 


67553 


0.23 


0.35 


0.26 


8460 




0.30 


-0.01 


907 


0.26 


0.31 


0.37 


Sun:our linelist 




2.25 


2.34 


0.52 



Sun:Grevesse c& Sauval (1998) 
224 2J3 OST 



UVES 



ID 


[Na/Fe] 


[Y/Fe] 


[Zr/Fe] 


19925 


0.51 


0.43 


0.24 


20766 


0.53 


0.35 


0.33 


21191 


0.51 


0.38 


0.25 


21728 


0.37 


0.41 


0.35 


22089 


0.50 


0.34 


0.31 


24590 


0.30 


0.58 


0.32 


25709 


0.34 


0.46 


0.31 


26794 


0.36 


0.57 


0.44 


27448 


0.11 


0.34 


0.26 


28103 


0.17 


0.40 


0.09 


28356 


0.37 


0.49 


0.34 


28797 


0.44 


0.45 


0.43 


28847 


0.08 


0.49 


0.29 


28977 


0.40 


0.64 


0.46 


29027 


0.02 


0.47 


0.23 


29065 


0.17 


0.36 


0.24 


29222 


0.24 


0.57 


0.28 


29272 


0.05 


0.47 


0.07 


29282 


0.42 


0.44 


0.35 


29545 


-0.02 


0.43 


0.25 


29598 


0.40 


0.53 


0.25 


29848 


0.09 


0.53 


0.17 


30209 


-0.05 


0.43 


0.13 


Sun:our linelist 




6.32 


2.25 


2.56 


Sun:Grevesse & Sauval_(1998) 




6.33 


2.24 


2.60 
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Table 5. Estimated errors on abundances due to errors on atmospheric parameters and to spectral noise compared with the observed 
errors 



ID 


ATeff=+50 K 


Alog(g)=+0.10 


A[Fe/H]=+0.05 


Av,=+0.10km/s 


S/N 


CTtot 


T obs 


A(loge(Li)) 


+0.05 


0.00 


0.00 


0.00 


0.07 


0.09 


0.10 


A([C/Fe]) 


-0.02 


0.01 


0.03 


0.03 


0.05 


0.07 


0.07 


A([N/Fe]) 


0.05 


0.01 


0.02 


0.03 


0.05 


0.08 


0.09 


A([0/Fe]) 


-0.04 


0.08 


0.01 


0.03 


0.06 


0.11 


0.10 


A([Na/Fe]) 


-0.02 


0.00 


0.00 


0.02 


0.08 


0.08 


0.08 


A([Mg/Fe]) 


-0.02 


-0.02 


0.01 


0.00 


0.08 


0.09 


0.07 


A([A1/Fe]) 


-0.02 


0.01 


0.00 


0.03 


0.07 


0.08 


0.09 


A([Si/Fe]) 


-0.04 


0.02 


0.01 


0.02 


0.06 


0.08 


0.06 


A / r /"^ rr~' ~i \ 

A([Ca/Fe]) 


0.00 


0.00 


0.00 


0.01 


0.04 


0.04 


0.04 


A([Ti/Fe]) 


0.02 


0.00 


0.00 


0.01 


0.03 


0.04 


0.04 


A([Cr/Fe]) 


0.00 


0.00 


0.00 


0.02 


0.07 


0.07 


0.08 


A([Fe/H]) 


+0.05 


-0.01 


-0.01 


-0.03 


0.01 


0.06 


0.05 


A([Ni/Fe]) 


-0.01 


0.01 


0.00 


0.01 


0.03 


0.03 


0.03 


A([Y/Fe]) 


-0.03 


0.06 


0.01 


-0.08 


0.06 


0.12 


0.14 


A([Ba/Fe]) 


-0.02 


0.03 


0.03 


-0.02 


0.04 


0.06 


0.04 


A([Eu/Fe]) 


-0.05 


0.04 


0.01 


0.03 


0.07 


0.10 


0.10 


A('2C/'3C) 














3 


3 


4 



Table 6. Equivalent Widths. Table|6]is published in its entirety in the electronic edition of the Journal. 



Wavelength(A) 


Element 


E.P.(eV) 


log(gf) 


#28590 


5711.083 


12.0 


4.34 


-1.67 


86.4 


6696.014 


13.0 


3.14 


-1.56 


34.7 


6698.663 


13.0 


3.14 


-1.83 


14.3 


5645.603 


14.0 


4.93 


-2.12 


23.6 


5690.419 


14.0 


4.93 


-1.84 


35.0 


5793.066 


14.0 


4.93 


-2.02 


28.4 


6125.014 


14.0 


5.61 


-1.58 


14.3 


6145.010 


14.0 


5.61 


-1.45 


20.7 


6244.465 


14.0 


5.62 


-1.34 


23.9 


6161.287 


20.0 


2.52 


-1.29 


54.3 


6162.170 


20.0 


1.90 


0.46 


170.9 


6166.429 


20.0 


2.52 


-1.14 


53.2 


6126.214 


22.0 


1.07 


-1.36 


26.8 


6258.098 


22.0 


1.44 


-0.34 


49.7 


6261.094 


22.0 


1.43 


-0.44 


42.4 
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